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Abstract 

We derive an exact functional renormalization group equation for the projectable version of 
Hof ava-Lifshitz gravity. The flow equation encodes the gravitational degrees of freedom in terms 
of the lapse function, shift vector and spatial metric and is manifestly invariant under back- 
ground foliation-preserving diffeomorphisms. Its relation to similar flow equations for gravity in 
the metric formalism is discussed in detail, and we argue that the space of action functionals, 
invariant under the full diffeomorphism group, forms a subspace of the latter invariant under 
renormalization group transformations. As a first application we study the RG flow of the New- 
ton constant and the cosmological constant in the ADM formalism. In particular we show that 
the non-Gaussian fixed point found in the metric formulation is qualitatively unaffected by the 
change of variables and persists also for Lorentzian signature metrics. 



1 Introduction 



The construction of a consistent and predictive quantum theory for gravity constitutes one of the 
major chahenges for theoretical high energy physics to date. Within the realm of quantum field 
theory there are currently two proposals that receive a lot of attention. Quantum Einstein Gravity 
(QEG) [1, 2, 3, 4] is based on Weinberg's Asymptotic Safety scenario [5, 6] and suggests that 
gravity could be a non-perturbatively renormalizable quantum field theory. In this case the UV- 
completion of the theory is provided by a non-Gaussian fixed point (NGFP) of the renormalization 
group (RG) flow. For RG trajectories that are captured by the NGFP in the UV, this construction 
ensures that physical quantities are free from unphysical UV divergences. Provided that the UV 
fixed point comes with a finite number of relevant directions, this construction is predictive. 
Starting from the seminal work [7] there is now a solid body of evidence supporting the existence 
and predictivity of this NGFP [8, 9, 10, 11, 12, 13, 14, 15, 16, 17]. The more recent proposal 
advocated by Hofava [18] suggests that the UV completion of gravity is an action which is 
anisotropic in the space and time coordinates. Higher powers of the spatial momenta entering 
into the gravitational propagators could render this construction perturbatively renormalizable 
[19, 20] but lead to difficulties when trying to restore Lorentz invariance in the low energy regime, 
see [21, 22, 23] for selected reviews and further references. 

A crucial ingredient in investigating both scenarios is the RG flow of the theory. For prac- 
tical computations, this flow can conveniently be encoded by the Wetterich equation [24] which 
captures the scale dependence of the effective average action F^ and can schematically be written 
as 



Here, F^ denotes the second functional derivative of the effective average action with respect 
to the fluctuation fields (p at fixed background (f) and Tr contains a sum over all fields of the 
theory and an integration over loop-momenta. Furthermore, TZk is a matrix-valued IR cutoff, 
which provides a A;-dependent mass term for fiuctuations with momenta < k"^. The interplay 
between the full regularized propagator ^F^^^ -|- TZk^ and kdkR-k ensures that the Tr receives 
contributions from a small p^-interval around A;^ ~ only, rendering the trace contribution finite. 

The flow equation (1.1) is deflned on the so-called theory space. This space contains all action 
functionals that can be build from a given fleld content and are compatible with the symmetries of 
the theory under consideration. Considering metric QEG, for example, the theory space T'^QEG 




(1.1) 
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contains all action functionals build from the spacetime metric g^j,, which are invariant under 
general coordinate transformations Diff(A^). 

An interesting open question at this stage is if the metric g^i, indeed provides the correct 
formulation for the gravitational degrees of freedom. Coupling gravity to fermionic matter degrees 
of freedom suggests that a first-order formalism based on the vielbein may be more fundamental. 
A first investigation of this scenario was initiated in [25] utilizing the Wetterich equation tailored 
to the fundamental fields of Einstein-Cartan gravity, i.e., the vielbein and the spin connection. 
Along a different line, the existence of time points at a preferred direction which in terms of 
Euclidean geometry may reflect itself in a foliation structure of spacetime. Such a structure is 
naturally captured by the ADM decomposition of the metric [26]. While all these formulations 
may be on-shell equivalent at the classical level, it is not clear if they describe the same quantum 
theory. In particular it is a priori unclear if the NGFP underlying Asymptotic Safety in the metric 
formulation also appears when using different fields for encoding gravity and, if so, weather these 
descriptions fall into the same universality class, in the sense that the universal critical exponents 
associated with the NGFPs actually coincide. 

In this paper, we address these issues by adapting the Wetterich equation (1.1) to the ADM 
decomposition, thereby implementing a foliation structure on the underlying quantum spacetime. 
In this case the metric g^i, is decomposed into a lapse function A^, a shift-vector Ni and a metric 
on spatial slices aij 



The lapse function and the shift vector essentially describe how the spatial slices are welded 
together, thus imprinting the (Euclidean) spacetime M with a preferred direction. This decom- 
position then naturally entails a spacetime structure that is topologically Ai = x Ti, where S 
are the d-dimensional leaves of the foliation. The direction singled out by the allows to Wick- 
rotate between Euclidean and Lorentzian signature, and will thus be referred to as (Euclidean) 
time direction. Notably, the ADM construction is very close to the geometric setting underlying 
the Monte-Carlo simulations of Causal Dynamical Triangulations (CDT) [27, 28, 29, 30, 31]. 

Depending on the precise field content and symmetry groups there are several theory spaces 
that can naturally be associated with the decomposition (1.2), see Table 1. Firstly, one can 
insist that the symmetry group acting on the ADM fields is the full Diff(A^) symmetry. In this 
case eq. (1.2) constitutes a non-linear field-redefinition of the gravitational degrees of freedom. 
The resulting theory space of foliated QEG, 7~^Q^*^^ is equivalent to the one of metric QEG, 



g^j,u '-^ {N , Ni, aij } . 
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theory 


theory space 


gravitational fields 


symmetry 


metric QEG 


-y-mQEG 




Difr(A^) 


foliated QEG 
projectable Hofava-Lifshitz 
non-projectable Hofava-Lifshitz 


■y-fQEG 
jT>HL 
-y-npHL 


N{t, x) , Ni{T, x) , aij{T, x) 
N{t) , Ni{T,x) , cry(r, x) 

N{t, x) , Ni{T, X) , (JijiT, x) 


Difr(A^) 
Diff(7W,E) 
Diff(M,S) 



Table 1: Definition of the gravitational theory spaces T that emerge within the ADM decompo- 
sition of the metric field (1.2). According to their precise field content and underlying symmetry 
group they constitute natural generalizations of the theory space underlying metric QEG, see 
(1.3). 

-y-mQEG^ j^g interaction monomials are constructed from the ADM fields and preserve Diff(A^). 
The foliation structure of Ai naturally defines a subgroup of Diff(Al), the foliation-preserving 
diffeomorphisms Diff(7W,S), eq. (3.8). Adopting this subgroup as the symmetry group of the 
theory space gives rise to the theory spaces underlying Hofava-Lifshitz (HL) gravity [18]. In 
its projectable version (pHL), the lapse function N{t) depends on time only. In terms of the 
field content, this can be understood as a partial gauge fixing eliminating the space dependence 
of the lapse function living on T^Q^G!^ The weaker symmetry requirements on 7"?^^ q^j^q-^ 
write additional interaction monomials for Lfc, which are invariant under Diff(A^,S) but break 
Diff(A^) invariance. Thus T^Q^*^ is embedded into the theory space of projectable HL theory 
j-pHL Finally, the non-projectable version of HL theory also includes a spatial dependence in 
the lapse function. Thus 7""?^^^ contains additional interaction monomials not present in 7^^^^ 
Based on these different symmetry requirements the theory spaces thus satisfy 



j-mQEG _ j-fQEG ^ -y-pHL ^ -y-npHL 



(1.3) 



As it will turn out, the natural symmetry of the Wetterich equation formulated in terms of 
the ADM fields are foliation-preserving diffeomorphisms. The key observation underlying this 
assessment is that Diff(7W) acts non-linearly on the ADM fields. Since it is a key requirement in 
the derivation of a Wetterich-type flow equation that its regulator is quadratic in the fluctuation 
fields, it is impossible to retain invariance with respect to a non-linearly realized symmetry. Thus 
our functional renormalization group equation (FRGE) is invariant under background foliation- 
preserving diffeomorphisms only, since Diff(A^,S) is the maximal subgroup that is realized by 
linear transformations. This restriction of the symmetry group has a direct consequence for the 
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theory space on which our FRGE is formulated: besides containing all interaction monomials that 
are invariant under the DifF(A^), the theory space also contains interactions which are invariant 
under foliation preserving diffeomorphisms and thus break the full diffeomorphism invariance 
explicitly. Thus our flow equation naturally encodes the RG flow on T'pHL i 

The relation (1.3) indicates that T'niQEG ^^^^ embedded into fP^^. An important result 
following from [7] is that the subspace X^QEg; jg actually closed under the RG flow. If the 
RG flow starts from an action functional preserving full diffeomorphism invariance, p^Q^*^ g 

k 

-y-niQEG^ integrating out quantum fluctuations will not generate interactions that violate Diff(A^) 
dynamically. This leaves the phenomenologically interesting scenario that T'^QEG ggj^-yes as an 
IR-attractor within T'pHL^ leading to a dynamical restoration of Lorentz symmetry at low energies. 

The foliation structure 7W = 5^ x S will actually lead to a FRGE which closely resembles 
the ones encountered for quantum field theories at finite temperature. In this interpretation, the 
system is seen as d-dimensional gravity coupled to a scalar N and vector Ni at finite temperature 
T proportional to the radius R of the S^. In this setting quantum fluctuations are split into 
two classes: spatial fluctuations of the d-dimensional system at ii = (zero temperature) and 
thermal fluctuations associated with the circle 5^. The latter are captured by the Matsubara 
modes of the system. The FRGE constructed in this paper integrates out both quantum and 
thermal fluctuations, essentially relying on the "imaginary time" formalism of finite temperature 
quantum field theory. This analogy allows to parallel the construction of Wilsonian flow equations 
for thermal fleld theories [32, 33, 34], adapting it to the gravitational system at hand. 

The rest of this work is organized as follows. We start by reviewing the classical ADM 
construction [26] in Sect. 2. The Wetterich equation capturing the RG flow on 7"?^^ jg constructed 
in Sect. 3. This flow equation respects foliation-preserving diffeomorphisms as a background 
symmetry, which severely restricts the interaction monomials that can be generated dynamically 
by the RG flow. As a flrst application we consider a truncation of the effective average action F^ 
which does not contain anisotropic couplings and derive the signature-dependent beta functions 
of the ADM-decomposed Einstein-Hilbert truncation in Sect. 4? Their properties and flxed point 
structure, which have partially been reported previously [35], are studied in in Sect. 5. We close 
with a summary and an outlook on possible applications in Sect. 6. The technical details entering 
the construction of the beta functions have been collected in two appendices. 

^As it will turn out, the off-shell setting of the FRGE implies that foliated QEG and projective Hofava-Lifshitz 

gravity are actually described by the same flow equation. 

^First results on the RG flow on fT^^ including Diff(A^)-breaking interactions will be reported elsewhere [36]. 
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2 Decomposing spacetime into space and time 

In order to make the present work self-contained and fix our notation, we start by reviewing the 
classical ADM decomposition of the metric field [26]. This construction plays an essential role 
when studying gravitational RG flows on spacetimes carrying a foliation structure or exhibiting 
anisotropic-scaling effects between space and time. 

The ADM construction starts from a D-dimensional Riemannian manifold Ai with metric . 
The spacetime signature can either be Euclidean, (+, +,...), or Lorentzian, (— , +,+,.. .) and we 
use the signature parameter e = sig7Q,^ = ±1 to distinguish the two cases. On this manifold we 
introduce a time function t{x) relating a real number, which we call time r, to every spacetime 
point X. This function equips A4 with a vector field SqT which, for Lorentzian signature, is taken 
to be timelike. The D-dimensional spacetime manifold can then be seen as a stack of spatial 
slices = {x '■ t(x) = Tj} with spatial dimension d = (D — 1). The hypersurfaces S^-- come 
with a normal vector Ua, which WG take as future-directed (^ri^dc^ r > 0) and normalized to unity 
{'^^^rLaUj^ = e). It can be related to the time function by = eNdaT where the Lapse function 
N acts normalization factor. 

On each spatial slice we introduce coordinates y*, i = 1,. . . ,d. On neighboring slices these 
coordinate systems will be related by the integral curves along daT. Explicitly we choose to 
be constant along such a curve. The Jacobian relating the coordinate systems x^ and (r, ?/*) is 
given by 

, ef=S . (2.1) 



dr 



_ dx" 



Note that n^ef = holds, since is normal to the hypersurfaces. The vector can be 
decomposed into its components tangential and perpendicular to the spatial slices as 

e = Nn'^ + N'ef . (2.2) 

Here, the d-dimensional shift vector A^* is purely spatial. 

For the coordinate one- forms the change of coordinates (2.1) implies 



dx" = edr + efdy' = Nn^dr + \^dy' + N'dr ) ef . (2.3) 
Consequently, the infinitesimal squared line element is given by 

ds"^ = japdx^dx'^ = eN^dT^ + dij idy' + N'dr) ( dy^ + N^dT] , (2.4) 
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where aij = ^^15 efe^ is the induced metric on the spatial shces and the component fields depend 
on the spacetime coordinates 

^fQEG = { jv(r,x^) , N,iT,x') , a,,{T,x') }. (2.5) 

From this expression we read off' the relation between the spacetime metric 7q,^ and the ADM 
fields 

Here the scalar products are with respect to the spatial metric aij. For completeness, we note 
that the decomposition of the determinant, appearing in the spacetime volume, is given as = 

3 The RG equation for foliated spacetimes 

In the sequel, we will adopt the viewpoint that the gravitational degrees of freedom are carried 
by the Lapse function N, the shift vector Ni and the spatial metric aij and derive the Wetterich 
equation for the component fields. We start with a discussion of symmetries and possible gauge 
fixings in subsection 3.1 before deriving the actual fiow equation in subsection 3.2. This construc- 
tion will lead to a FRGE which intrinsically implements the foliation structure of spacetime and 
captures the RG flow of the projectable case of Hofava-Lifshitz gravity [18].^ 

3.1 Symmetries and gauge fixing 

Since symmetries play a crucial role when constructing a FRGE, we start with a systematic 
discussion of the diffeomorphism symmetry in the ADM framework. 

3.1.1 Classical gravity in the ADM formalism 

Under a general coordinate transformation Diff(A^) the spacetime metric transforms according 

to 7a/3 7q,/3 + 6-/al3 with 

h^lu = C-vl^lv ■ (3.1) 



■^As discussed in the introduction, the theory space underlying the non-projectable version of Hoi^ava-Lifshitz 
gravity admits an additional class of interaction invariants constructed from the new building block ai — diN/N. 
Since in this case there are additional subtleties concerning the closure of the constraint algebra [37, 38] it will not 
be discussed further at this point. 
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Here denotes the Lie derivative of 7 with respect to the D-dimensional vector v°^{x^). The 
decomposition (2.2) then ahows to write v°' in terms of its time and spatial parts f{T,x) = f^Va 
and C{t, x) = e^w" 

^" = (/(r,x),C(T,x)). (3.2) 

Combining this spht and the ADM decomposition (2.6) ahows us to determine the transformation 
behavior of the component fields under Diff(A^) 

6N = drifN) + C'^dkN - NN'd^f , 

6Ni = dr{Nif) + C'^dkNi + iVfca.C' + ^fc^a,C' + NkN^'dif + eN^d^f , (3.3) 
daij = fdr^ij + C^dkdij + djkdiC^ + ^ikdjC^ + Njdif + Nidjf . 

For completeness, we note 

6N' = dr{N'f) + C^djN' - N^djC + drC - N'N^d^f + eN^a'^djf . (3.4) 

One observes that, while Diff(A4) acts linearly on the metric 7^,^, the non-linearity of the ADM 
decomposition (2.6) leads to a non-linear transformation law for the component fields. 
The gauge freedom can be exploited to adopt the proper-time gauge [39] 

N{t,x) = 1, N\t,x) = 0. (3.5) 

This gauge choice still leaves the freedom to choose coordinate transformations which satisfy 
drf = and drC + ^'^^'^dkf = 0. Note that these equations encode the freedom to choose a 
coordinate system on the initial slice of the foliation. At the level of the path integral they are 
typically fixed by the corresponding boundary conditions, see [40, 41] for a more detailed discus- 
sion. Since we are not interested in surface effects, we will neglect this point in the subsequent 
discussion. 

3.1.2 Hofava-Lifshitz gravity: the projectable case 

Hofava-Lifshitz gravity [18] encodes the gravitational degrees of freedom in terms of the ADM 
fields. Considering the projectable version of the theory, the key difference is that here only the 
metric on the spatial slices and the shift vector are spacetime dependent fields, while the lapse 
function N[t) depends on time only and is constant along S 

;^PHL = |^(^)^ ^^(^^^)^ ~..^^^^)|, (3g) 
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Moreover, the symmetry group is restricted to foliation-preserving diffeomorphisms Diff(A4,S). 
In this case the vector appearing in the Lie derivative (3.1) is restricted to the form 

T;°(r,x) = (/(r),C(r,x)) . (3.7) 

The variations of the component fields under foliation-preserving diffeomorphisms are those 
of (3.3) restricted by the fact that / and N are space independent 

SN = drUN) , 

5Ni = drifNi) + C^dkNi + NkdiC"" + aikdrC"" , (3-8) 

while 

SN^ = drifN') + drC - N^dkC + ('"dkN' ■ (3.9) 

In contrast to the full diffeomorphisms (3.3) the foliation-preserving diffeomorphisms act lin- 
early on the component fields. This will lead to considerable simplifications when applying the 
background-field method later on. 

Again, we can use the gauge freedom to adopt the proper-time gauge (3.5). This fixes the 
Diff(A^,S) up to the residual transformations drf = and drC'' = 0, i.e., the choice of co- 
ordinate system on the initial slice. Thus, upon gauge-fixing, the field content of projectable 
Hofava-Lifshitz gravity and foliated QEG are identical. In this sense the transition from the field 
content (2.5) and gauge-symmetries (3.2) of the classical ADM formalism to (3.6) and (3.7) can 
be understood as a partial gauge-fixing of the former. Owed to their identical off-shell field con- 
tent, the Wetterich equation for the ADM formulation of classical gravity T^Q^^; ^^^^^ projectable 
Hofava-Lifshitz gravity fP^^ will look identical. We stress, however, that the two theories are dif- 
ferent in the sense that the latter allows a larger class of admissible interaction functionals, since 
the requirement of invariance under foliation-preserving diffeomorphisms is less restrictive then 
demanding invariance under Diff(A^). At this stage we proceed by constructing the Wetterich 
equation on TT^l^ which retains Diff(Al,S) as a background symmetry. 

3.1.3 Gauge fixing in the background field formalism 

The consistent quantization of the theory requires gauge fixing the symmetries (3.8) in order to 
restrict the integration in the path integral to physically inequivalent configurations. For our 
purpose it is most convenient to implement this gauge fixing via the background field method 
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[42], see [1] for a detailed review. This construction ensures that the resulting effective action 
contains interaction monomials that are invariant under (3.8) only. 

When implementing the background field method the quantum fields (3.6) are split into 
fixed but arbitrary background fields x = , ^i, (^ij} and fluctuations around this background 

N = N + N, Ni = Ni + Ni, aij = aij + aij. (3.10) 

Note that the fluctuation flelds are not assumed to be small in any sense. No expansion in powers 
of X is implied in this split. 

The central element of the background fleld method is that the symmetry transformations 
X '"^ X + "^X can be distributed between the background and fluctuation fields in different ways. 
Quantum gauge transformations leave the background invariant and attribute (3.8) completely 
to the fluctuation field, 

x^x, x^x + '^q(x + x)- (3.11) 

Here the variations 6q are given by (3.8) with x substituted by x + X- It is this set of symmetries 
that have to be gauge-fixed. 

Background gauge transformations on the other hand act on both the background and fiuc- 
tuation fields 

X^X + ^BX, X^X + f^BX- (3.12) 

This transformation plays the role of an auxiliary symmetry, which is retained by all the terms 
entering into the path integral. Its purpose is to ensure that the effective action contains interac- 
tion monomials that are invariant under Diff(A^, S) only. It thereby suffices that the background 
transformations become identical to the ones of the quantum field once the fiuctuations are set 
to zero. This leaves some freedom in the choice of 5b- Explicitly, we adopt 

5bN = dr ifN) , 5BN^ = dr {f Ni) + aijOrC + , SB^7^j =fdra,j + C^{aij) (3.13) 

for the background flelds while the fluctuations are taken to transform as 

SbN = dr (/iv) , SbN., = dr (/iV,) + C^{Ni) , 5b<7.j =fdra^j + C(;{aij) . (3.14) 

Here £^ denotes the Lie derivative on the spatial slices which contains only spatial derivatives. 
Note that here the drC-^erms have all been incorporated in the transformation of the background 
field. While this implies that the background fields in (3.13) do not transform as the corre- 
sponding tensors, this still constitutes a choice of the background split since it reduces to the 
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transformations (3.8) in the limit of vanishing fluctuation fields. Combining these transformation 
laws shows that 6-bX is given by (3.8). Here it is essential that Diff(A^,S) acts linearly on x- 

In the next step, we construct a gauge fixing term that implements the proper-time gauge 
in the background field formalism. By definition, this term has to break the quantum gauge 
transformation (3.11) while being invariant under the background gauge transformation (3.12). 
A straightforward computation then establishes that 

iV2 1 



= \V^J dT(fxN^/d 



+ — (Jij 



ai iV2 as 



(3.15) 



indeed satisfied these requirements. In the Landau-limit, — )• 0, 05 — )• 0, S^^ becomes a delta- 
distribution which eliminates the fluctuations of the lapse and shift vector, implementing the 
background proper-time gauge 

N ■ 

i^lapse = ]^ = , i^shift ^ ^ = • (3.16) 

The ghost action, exponentiating the Faddeev- Popov determinant arising from (3.15) can then 
be found in the standard way 



y/e J drd'^x \fo 



(3.17) 



The vector ghosts u)i{T, x),uj^{t, x) are functions of the spacetime, while the scalar ghosts Cj{t),uj{t) 
depend on time only. The background lapse function N{t) has been distributed in such a way 
that the ghost action is invariant under background gauge transformations with the ghosts trans- 
forming as scalars and vectors, respectively, 

6bU} =f drUJ , 6BU!i = f drOi + £^(wj) , 5buj^ = f drUJ^ + Cc^{uj^) . (3.18) 

3.2 The Functional RG Equation 

After discussing the symmetries of our gravitational theory, we are now in the position to derive 
the Wetterich equation encoding the RG flow on T'pHL^ The construction follows the standard 
derivation [24, 7], see, e.g., [43, 44] for pedagogical reviews in the context of gravity. 

3.2.1 Defining the eflPective average action 

Our starting point is a generic action S^^^^[N ^ Ni, dij] build from the multiplet (3.6) and invariant 
under foliation-preserving diffeomorphisms (3.8). We then consider the sca/e- dependent generating 
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functional Wk for the connected Green functions 



exp{Wk[J; x]} = jvii exp \^-S^''-'' - S^^ - - A^^ - 5'^°"''™] . (3.19) 

Here S^^ and 5"^^ are given by (3.15) and (3.17), respectively and the measure T>^ consists of 
the integration over the gravitational fluctuations DNDNi'Daij and the ghost contributions. The 
source term 



^source = _ / rfrd^'xiV 



t ^ - • • - - • • 



(3.20) 



is invariant with respect to background-Diff(7W, S) if the sources transform covariantly with 
respect to time-reparametrizations / and via the Lie derivative with respect to spatial diffeo- 
morphisms Q. Here indices are raised and lowered with the spatial background metric. For 
convenience, we absorb the extra powers of N appearing in the lapse and shift terms into the 
sources, setting t = tN^'^ and = t*iV~^. The sources can then be collectively written as 

J = {t,t\f^,^,^,^i,^'}. (3.21) 

The essential piece in eq. (3.19) is the IR cutoff for the gravitational multiplet and the ghosts 

(3.22) 

+ yfe drd^xNVd^ {Cd,idi) TZfix] {uj,ujjf ■ 



The cutoff operators TZ^^^ [x] and TZ^ [x] serve the purpose of discriminating between the high- 
momentum and low-momentum fluctuations. Following [7], we use the eigenvalues of the Lapla- 
cian constructed from the background multiplet x to discriminate these modes. Eigenmodes of 
—D^ with eigenvalues ^ k"^ are integrated out without suppression whereas modes with small 
eigenvalues <^ k"^ are suppressed by a momentum dependent mass term. Generally, the TZk 
have the structure 

IZM = ZkRk = Zkk^R^^H-D^/k^) . (3.23) 

Here Zj, is a background-field dependent matrix that ensures the invariance of (3.22) with respect 
to the background-gauge transformations. The dimensionless shape function R^^") interpolates 
between i2^^)(0) = 1 and Imiw^oo R^^\w) = 0. Convenient choices are, e.g., the exponential 
cutoff R^^\w) = w[eyiY>{w) — 1]~^ or the optimized cutoff (A. 12). 
In the next step, we construct the /c-dependent classical fields 

X = {h, hi, hij , c, c, Ci, d} (3.24) 
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as the expectation value of the fluctuation fields. These are easily found as variations of Wk with 
respect to the sources 

y/eN^/a oJa 

The extra factors N included in (3.21) ensure that the classical lapse and shift transform as the 
corresponding quantum fields 

dBh=dr{fh), SBhi = dr{fhi)+C^{hi). (3.26) 

As usual, we assume that one can invert the relations (3.25) and solve for the sources as 
functionals of the classical fields and, parametrically, of the background. The Legendre transform 
Tfc of Wk reads 

f x] = Vej dT(fxNV^ [th + fhi + f^hij +^c + dc + ^ic' + i^'ci] - Wk . (3.27) 

The effective average action Tk is then obtained from Tk by subtracting the cutoff action AkS 
with the classical fields inserted 

Tfc [x;x] = f [x; x] - A,.5[x; x] ■ (3.28) 

3.2.2 The Wetterich equation 

The derivation of the Wetterich equation encoding the /c-dependence of starts with the con- 
nected two-point function 

iGkr\x,T,-y,Ty) = ^ = ^'Y/ -, (3.29) 

~ (2) 

which is matrix-valued in field space. Correspondingly, the Hessian of (3.27) is given by 

^''^ab- e iV2 y^a{x, T,) a{y, Ty) Sx''{x, rj<5x^(y, Ty) ' 

where the index [b] takes the value zero for commuting and one for anti-commuting fields Xfti 
respectively. The Legendre transform (3.27) implies that (3.29) and (3.30) are each others inverse 
in a functional sense 

I d'ydTyNV^iGkr{x,T,;y,Ty) (ff r,; z, r,) = . (3.31) 

Taking the derivative of (3.19) with respect to the RG "time" t = ln(k) then gives 

1 



-dtWk = -TT 



{r^XldtUkU . (3.32) 
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Here {x°'^x^) denotes the two-point correlator of the fluctuation fields (including the ghost fields) 
and Tr includes an integration over loop momenta and a sum over internal indices. The two-point 
correlator is related to the connected two-point function {Gkjab ^-nd the classical fields, (3.25), 
via 

{r ® X") = {Gkr' + = (f + ■ (3.33) 

\ / ab 

Expressing the flow (3.32) in terms of the effective average action (3.28) then yields the desired 
functional renormalization group equation for Hofava-Lifshitz gravity 



dt^k[x;x] = ^STr 



(3.34) 



Here the Hessian ' denotes the second derivative of F^fx; x] with respect to the fields x- Both 

(2) 

and TZk are matrix valued in field space and the supertrace includes an integral over loop 
momenta and a sum over field space. The FRGE (3.34) constitutes the central result of this 
section. 

We close this subsection by discussing the limits of the RG flow implied by (3.34). Since 
the TZk's vanish for k = 0, the limit A; — )• of rfc[x;x] brings us back to the standard effective- 
action functional which still depends on two sets of fields x ^^d x- The ordinary effective action 
r[x] with one argument is obtained from this functional by setting the expectation value of the 
fluctuation fields x = [42] 

r[x] = limrfc[x = 0;x]. (3.35) 

Besides the FRGE (3.34) the effective average action also satisfies an exact integro-differential 
equation, which can be used to find the A; — )• oo limit of the average action [7]: 

rfc^oo[x;x] = s[x + x] + 'S'gf[x;x] + 'S'gh[x;x] • (3-36) 

Intuitively, this limit can be understood from the observation that for A; — )• oo all quantum 
fluctuation in the path integral are suppressed by an infinite mass term. Thus, in this limit no 
fluctuations are integrated out and T^^oo agrees with the microscopic action S supplemented by 
the gauge fixing and ghost actions, also see [45] for more details. 

3.2.3 Symmetries preserved by the RG flow 

The main advantage in the use of the background field method is that the functional F^ and, as 
a result, also F^ is invariant under background foliation-preserving diffeomorphisms, when all its 

14 



arguments transform according to their transformation rules 

^k[x + SBx;x + SBx] = '^k[x;x]- (3-37) 

Note that here, contrary to the quantum-gauge transformation (3.11), also the background fields 
transform according to their corresponding symmetries. The invariance (3.37) is a consequence 
of 

Wk[J + 6BJ;x + SBx] = Wk[J;x] (3.38) 

which in turn follows from (3.32) if one uses the invariance of S^"^^^, S^^, S^^ and A^S under 
the background transformations (3.12) and (3.18). At this point, we assume that the functional 
measure in (3.32) is invariant under Diff(Al,S). 

The background-gauge invariance of Ffc, expressed in eq. (3.37), is of enormous practical 
importance. It implies that a RG flow, starting from a background foliation-preserving diffeo- 
morphism invariant at a scale k, will not generate interactions that violate this background 
symmetry dynamically. Nevertheless, even if the initial action is simple, the RG flow will generate 
all sorts of local and non-local terms in which are consistent with the symmetries. 

We close this section with the following important remark. Owed to the requirement that the 
cutoff AfcS" must be quadratic in the fluctuation fields, it generically seems impossible to construct 
a Wetterich-type FRGE which employs a linear background split and preserves a non-linear 
symmetry as a background symmetry. In principle it should be possible, however, to construct 
a FRGE for the gravitational multiplet (2.5) where the full Diff(A^) is realized as a background 
symmetry. Since the corresponding symmetry transformations (3.3) are non-linear, we expect that 
this construction will either require a generalization of the background gauge fixing procedure 
to non-linear symmetries [46] or a non-linear background split [1]. From the flow equation for 
gravity in the metric formalism [7], we expect that this will lead to a functional renormalization 
group equation on a subspace T^Q^^; ^ -y-pHL .^^j^j^j-^ jg closed under RG transformations. While 
constructing the corresponding flow equation is certainly interesting we will not embark on this 
construction at this point. 

4 Signature effects in the gravitational RG flow 

In principle the Wetterich equation (3.34) constitutes an exact RG equation on T'pHL -which is, 
however, notoriously difficult to solve. A common technique to find approximate solutions of 
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the equation which do not rely on an expansion in a small parameter consists in truncating 
the effective average action by restricting T/. to a finite set of running coupling constants. In 
this section we will derive the beta functions in the simplest gravitational setting, the foliated 
Einstein-Hilbert truncation. This setup allows the direct comparison of the approximate RG flows 
on T'l^QEG ^j^^ -y-fQEG^ Owed to the foliated background, our derivation uncovers many structures 
that are well-known in the context of quantum field theory at finite temperature [32, 33, 34]. 

4.1 The ADM-decomposed Einstein-Hilbert truncation 

The ansatz for corresponding to the ADM-decomposed Einstein-Hilbert truncation is of the 
general form^ 

Tkix; x] « rr"[x + x] + s^'[x; x] + ^^^[x; x] (4.i) 

where we have approximated the gauge-fixing and ghost part of the effective average action 
by their classical expressions (3.15) and (3.17). For the gravitational part T^^^ we adopt the 
Einstein-Hilbert action in D = d+ 1 spacetime dimensions. Expressing 

= 16^ / ^""^^ + ^^'^ (^-^^ 
in terms of the gravitational multiplet (3.6) the corresponding action reads 

= I drd'^x [e-i {KijK'i - K^) - ^''^R + 2Afc] . (4.3) 

Here the {N,Ni, aij} are the classical counterparts of the full quantum fields x 

N = N + h, Ni = Ni + h,, a^j = aij + hij . (4.4) 

Moreover 

Kij = ^ {dra^, - DiN, - DjNi) (4.5) 

is the extrinsic curvature, K = a^^Kij, and Di and ^^^R denote the covariant derivative and 
the intrinsic curvature of the d-dimensional spatial slices constructed from aij, respectively. The 
functional form (4.3) can easily be established by substituting the ADM decomposition (2.6) 
into the standard Einstein-Hilbert action written in terms of the covariant metric 7q/3. As in the 
standard Einstein-Hilbert case (4.3) contains two scale-dependent coupling constants, the Newton 
constant Gk and the cosmological constant A^. 



^In the terminology of [47] this corresponds to a single-metric computation. Results for metric flows that take 
into account the effect of higher order terms in the fluctuation fields have recently been reported in [48, 49]. 
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The action (4.3) is, by construction, invariant under the full diffeomorphism group Diff(A^). 
This symmetry fixes the relative coefficients between the extrinsic and intrinsic curvature terms to 
the form (4.3). When restricting the symmetry group to the foliation-preserving diffeomorphisms 
(3.8), the term, the K^^ Kij term and the R term are invariant on their own. Thus the 
anisotropic case of Hof ava-Lifshitz gravity allows to introduce further coupling constants in front 
of the terms. The study of their RG flow will be subject to a future publication [36]. 



4.2 Constructing the functional traces 

Substituting the ansatz (4.1) into the FRGE (3.34), we observe from its left hand side that the 
scale dependence of and can be read off from the coefficients of the spacetime volume and 
the curvature on the spatial slices ^'^^R. Thus it is sufficient to project the traces appearing on 
the right hand side onto these two curvature invariants. We stress that when computing the beta 
functions of the theory the geometric quantities merely act as a bookkeeping devices. No physical 
meaning should be attached to them. 

The first step in evaluating the functional traces arising from (4.1) constitutes in computing 

(2) 

the Hessian [XjX] with respect to the fluctuation fields. This calculation can be simplified in 
a number of ways. Firstly, we adopt Landau gauge by sending — t- 0, — t- 0. In this limit the 
gauge-fixing term (3.15) is converted to a 5-function which freezes the fluctuation fields = 0, 
Ni = 0. Thus N and Ni decouple and do not contribute to the traces on the right hand side of 
the flow equation. Thus, 



dfT 



'ii k 



iTr 

2 



e2pgrav 

6hij5hki ^ 



Tr 



M+nf 



dtTZ 



(4.6) 



where the Hessian Ai in the ghost trace is given by 

M = diag 



(4.7) 



6c6c 6c^6c^ 

The second simplification originates from the observation that we can set the fluctuation fields 
X = once the second variations are computed. Thus we can use the background covariance of 
the construction to carry out the computation in a specific background geometry. The choice of 
background has to be general enough to distinguish the two interaction monomials carrying the 
information of the RG flow. A natural choice uses 



(4.8) 
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i.e., the direct product of a "time" -circle with periodicity T and a time-independent maximally 
symmetric sphere S'^. In terms of the background multiplet, this choice implies 

N{r) = l, N,{t,x) = 0, aij{T,x) = aij{x)\sci. (4.9) 

Owed to the maximal symmetry of S"^ the background curvatures on the spatial slices satisfy 



0. 



Ri 



R 



(4.10) 



^^J'^' - did-1) 

In order to lighten our notation, we will drop the prefix (d) of the intrinsic curvature from now 



on. 



We then proceed to compute the Hessian of Tf^^^ around the background (4.8). Inspecting 



(4.3), there are four distinguished interaction monomials entering into T^^^ 

h = jdTd'^xN^, l2 = j drd'^xN ^ R, 

I3 = / drd'^x N Vct KijK^i , 1^ = J drd'^x N . 



(4.11) 



Using the background (4.9), (4.10) to simplify the result, the second variations of these monomials 
read 



52 = J drd'^xV^i-^ h^'^ hjj + ^h^} , 



drd'^xVcri h 



(d-2)(d-l) A (d-2)(d-4) p 



h - i/l'^- 



^ + d(d^l) ^ 



^2/3 = - i y drd^xV^ [hi dl h^'^ + ^hd^,h], 
Here we decomposed the fluctuations hij into their traceless and trace part, 



(4.12) 



h 



T 1 - 

hij + -a-ij h , 



a^^hl = 



(4.13) 



and A = —a^^DiDj is the background Laplacian built from the spatial metric. Here and in 
the following the h denotes the trace-part of the classical fluctuation field hij and should not 
be confused with the classical expectation value of the lapse function, (3.25), which has been 
eliminated by the proper-time gauge choice (3.16). 

A technical difficulty arises from the fact that in (4.12) not all covariant derivatives combine 
into Laplace operators. In particular the second line in 6'^l2 leads to a so-called non-minimal 
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operator structure in r^^\^ Following [50, 14], we bypass this problem by further decomposing 
the traceless part of the metric fluctuations hjj into its irreducible representations on the sphere 
via the transverse-traceless (TT) decomposition [51] 

= h]^ + 2D^, [A - ii?] -^/^e,) + [A A- - 3^^.] [A(A - ^,R)] . (4.14) 

Here the component fields are subject to the differential constraints 

a'^hj^ = 0, &hj^ = 0, &Ci = 0. (4.15) 

Moreover, we have chosen the normalization of the transverse vector and the scalar ? such that the 
TT decomposition (4.14) substituted into the Gaussian integral f "Dhat exp[— J d'^xy/Bhijh^^] does 
not give rise to operator- valued determinants. Thus the normalized TT decomposition (4.14) does 
not require the introduction of auxiliary fields exponentiating the Jacobi determinants naturally 
appearing in this type of field redefinitions as, e.g., implemented in [52]. 

Substituting the TT decomposition (4.14) yields the second variation in terms of the compo- 
nent fields 

(4.16) 



1/2 



b^U = -\ \ drd'^xy/^h dl h . 

Here the coefficients Cj are functions of the dimension of the spatial slices d 

_ £-3d + 4 _ {d-2){d-l) 

d{d-l) ' = d^ • ^^-^^^ 

In order to derive (4.16) from (4.12) it is useful to first carry out the TT-decomposition without 
including the normalization factors for the transverse vector and scalar c and subsequently notic- 
ing that the normalization factors appear as factors in the hessians. Moreover, we used that the 
background metric is r-independent, such that all background quantities (including the spatial 
Laplacians) commute with the r derivatives, which leads to further drastic simplifications. 



^This feature was already expected based on earlier computations of the beta functions in the metric formulation 
[7, 50, 53], where a special gauge choice is required in order to remove these non-minimal terms. 
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Based on (4.16) it is now straightforward to write the part of r|'^^ quadratic in the fluctuation 



fields 



rr K-] = rr [^..] + o{h) + i ^^rf- a.,] . 



(4.18) 



It is convenient to list the contributions appearing in 6'^T^^''[hij; aij] according to the components 
of the fluctuation fields contained 



r2pgrav 




levrGfe 
IGttG, 

rzpgrav 

levrGfc 

i;2-pgrav \/^ 

^ ^f'^ "levrGfe 

^2pgrav 



C7o(A-i?)-^52- 



dTd'^xy/a\h 



2{d-l) 



(4.19) 



h, 



drd^'-xV^Coh 



1/2 



These variations are complemented by the analogous variations appearing in the ghost sector. 
Setting the background ghost fields to zero, which suffices for the present computation, (3.17) 
leads to second variations that are diagonal in the ghost fluctuations. Explicitly, they read 



(J^rg^^ = y drd'^xy/^cdrc, 
\fe / dTd^x\f^'S'drCi . 



(4.20) 



Based on the variations (4.19) and (4.20) we are now in a position to specify the IR regulators 
T^fc. As a key requirement Tt^ must be positive definite for all values of /c, in order to act as a 
/c-dependent mass term. In general, this will not be the case if A is a Laplacian constructed from 
a Lorentz signature metric. Since one of our interests is in the signature effects of the spacetime 
on the gravitational RG flow we choose a purely spatial regulator, i.e., 7^^ = 7^fc(A) will be 
a function of the Laplacian on the spatial slice only. In this case A is positive semi-definite. 
Moreover, this choice is sufficient to regularize the operator traces appearing in the flow equation 
[33]. 6 

In practice we will implement a regulator of Type I, in the nomenclature of [12], dressing the 
spatial Laplacians with a fc-dependent mass term according to 



A^ A + A:2 Rf{^lk^) 



k ■ 



(4.21) 



^Covariant cutoffs have recently been proposed in [34]. Since their implementation in the Lorentzian setting 
is slightly involved we will resort to the spatial regulators at this stage. It will be interesting to carry out the 
expansion around the mass poles of the Lorentzian propagators and compare to the results obtained here. 
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Here k"^ are purely spatial momenta and i?*-^) is the profile function introduced in (3.23). The 
rule (4.21) fixes the matrices Z}^ and we obtain 



T^TT 



1 



327rGfc 



n, 



327rGfc 
1 



(4.22) 



/it 



■Co 



p2 



1/2 



d-1 



RA 



1/2 



327rG/fc 

The quadratic variations of the transverse vector and the ghosts do not contain spatial Lapla- 
cians. Thus the corresponding regulators TZ^ constructed via (4.21) vanish. As a consequence, 
they do not source the RG flow of the Newton constant and the cosmological constant in the 
present truncation. 

We now have all the ingredients to write down the explicit form of the traces contributing to 
the flow of the truncation ansatz (4.3). Here it is convenient to separate the contributions of the 
transverse-traceless and scalar sectors according to 



dtTk = Ttt + Tsa + Ts,2 



(4.23) 



Explicitly, these traces are given by 



Ttt = ^ (327rGfc) TV [(Pfc - + C2tR - 2Afc) ' dtTlTT 



(4.24) 



and 



7;,i=i(327rGfc)Tr 



7;,2=i(327rGfc)Tr 



jf 2Co (P| - ^RPk 



(4.25) 



Here the cutoff functions are given in (4.22) and M is the determinant of the two-by-two matrix 
appearing in the scalar sector 

M =CoAk {Pk - 2Afc) - - 7^(3d - 4)A,92 



-%id-2)R( P, - m + ^R- ^A, 



(4.26) 

The traces and thereby capture the contribution from the diagonal and off-diagonal part 
of the scalar field matrix, respectively. 

For the computation of the RG flow, it is sufficient to keep track of the intrinsic curvature 
scalar up to the first order in R, since all higher powers of R are outside the truncation subspace. 
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This feature can be used to combine the two scalar traces by expanding the numerator in Ts^2 
around R = As a result, the two scalar contributions combine into 



% = %,i + Ts,2 ~ \ (327rGfc) Tr 



N 



2d 



+ 0{R^) (4.27) 



which considerably simplifies the evaluation of this sub-sector in the next subsection. 
4.3 Constructing the beta functions 

The projection of the traces (4.24) and (4.25) onto the truncation subspace conveniently exploits 
that the background has the structure x S'^. This structure allows to Fourier expand the 
fluctuations in the "time direction" according to 

oo T 

<A(r,x)= 0„(x)e2™-/^ ^ 4>^{x) = - dT<t>{T,x)e~^^'^^l^ (4.28) 

n=-oo -^0 

with the complex coefficients (pni^) obeying the reality constraint (f)n{x) = (/'(_„)(x)*. In this 
case the integration over the momenta in the r-direction, contained in Tr, is converted into a 
discrete sum over Matsubara modes. The remainder is given by tr which contains the operator 
trace restricted to the spatial slice together with a trace in field space. ^ 

oo 
n=— oo 

The r derivatives are then converted into Matsubara masses of the fiuctuations in time direction 

a^^-f^y. (430) 



T , 

The expansion of tr in terms of the intrinsic curvature can then be carried out by applying 
the standard heat-kernel expansion to the Laplace operator on the spatial slices. Applying the 
trace technology summarized in Appendix A.l it is then straightforward to find the two lowest 
terms in the spatial-curvature expansion. The result is most conveniently expressed in terms of 
the dimensionless analogs of Newton's constant, the cosmological constant and the Matsubara 
mass m 

2'jr 



^When performing a "dimensional reduction" of a higher-dimensional quantum field theory on a circle, the 
Matsubara modes are usually called Kaluza-Klein states. A string compactification typically restricts itself to the 
massless sector. In the language of field theory at finite temperature the latter corresponds to truncating the 
Matsubara sums considering the zero-mode only. 
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The dimensionality of Gk is the canonical one for a d + 1-dimensional spacetime. In addition, we 
introduce the anomalous dimension of Newtons constant 

7?^ = -dt \n{Tk ZNk) , Gk = Z~l Go . (4.32) 

Here we anticipated that the size of the time circle may depend on k, which would reflect itself 
through an anomalous running of the Matsubara mass m ^ k~^. 

The contribution of the transverse-traceless modes (4.24) then becomes 

T'tt = 1^ E / d'xV^k'^!^g]f,iw,^) + Rk~^ (i _,(^2t) - 4° (^2t)) } . (4.33) 

Here = \{d + ^){d — 2) arises from the trace over vector indices, 

W2T = e'^rn^n^ - 2Afc , (4.34) 

and the g-functions are defined in (A. 11). The scalar trace (4.27) can be evaluated along the 
same lines 

, d-2 r \ 2,-1/ N 1,0/ \ , (z 2 2 4(d-3) -> \ 2,0/ n] 

^2dx3 (^o) - Qd/sCw^o) + (^f m n - ) Q^/sC^^o) 

with 

^0 = (7^'"' - 2Afc) (im^n^ - ^A^) . (4.36) 

After working out the operator traces on the spatial slices the next step constitutes in per- 
forming the Matsubara sums for the "time- like" fluctuations. In order to be able to carry out 
the sums analytically, we specialize to the optimized cutoff (A. 12). In this case the thresh- 
old functions (A. 8) become particularly simple. Using the results derived in Appendix A. 2 the 
transverse-traceless part of the trace can be written as 

7-TT = ;|5^/Av/5fc-'{T;;M(K/l,-C,xT,y,)} (437) 

while the scalars contribute 

^ /■ ,d /^,d/c^l,0 1 i? c^l,0 ,d-2R\(. 4(d-3) , \ e2,0 ^.1,0 oo2,lll 

~ (ivrp/^ y " ^^^'^ l'^d/2 + 61^*^^/2-1 + 2dXn? [ y- d-2 ^d/2 + ^^d/2\ j ' 

(4.38) 

The functions Tj'''^ and Sf'"^ are defined in (A. 27) and their explicit forms are obtained in Appendix 
A.2. 
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The non- perturb at ive beta functions for the Newton constant and the cosmological constant, 
arising from the truncation (4.3), are encoded in the coefficients of the zeroth and first order 
terms of the intrinsic curvature. Taking the t derivative and setting X = X afterward, the l.h.s. 
of the fiow equation becomes 

a.r,|^^. = d^xV^{-Rd, (t) +20. . (4.39) 

Equating these curvature terms with the corresponding terms in (4.37) and (4.38) then gives rise 
to the beta functions for the dimensionfull Newton constant and cosmological constant. In terms 
of the dimensionless quantities (4.31) these read 

dt9k = Pg{gA]m) , dfXk = f3\{g,\;m) (4.40) 

where 

j3g = (d - 1 + tin) g , 



(4^)d/2 



/?, = ivN - 2)A + ( d2T T^f, + 



(4.41) 



The anomalous dimension of Newton's constant is given by 

mgBi{X) 
l + mgB2{\) 



VN = ""'""f],, (4.42) 



with 



+ ^((i-^a)^^;,-^^;,+3^^;,) 



(4^)d/2 



(4.43) 



The beta functions (4.40) constitute the final result of this section. Besides the dimensionless 
Newton constant and cosmological constant they parametrically depend on the Matsubara mass 
m, which encodes the size of the foliation in units of the IR-scale k. 



5 RG flow scenarios 

In this section we proceed by analyzing the RG flow emanating from the beta functions (4.40). 
Our prime interest thereby lies on the influence of signature effects on the existence of non-trivial 
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e 


A < -4Co 


-4Co < A < 1/2 


1/2 < A 


+1 


hyperbolic 


mixture 


trigonometric 


-1 


trigonometric 


mixture 


hyperbolic 



Table 2: The analytic structure of the beta functions (4.40) depends on the signature of spacetime 
e, the value of the cosmological constant, and the dimension of the spatial slice via the combination 
Co (4.17). We distinguish the cases where the threshold functions entering the beta functions are 
build from hyperbolic or trigonometric functions only or contain both types of terms. 

fixed points of the RG flow which could feature in a gravitational Asymptotic Safety scenario. 
Concretely, the general properties of the beta functions will be discussed in Sect. 5.1. Their fixed- 
point structure in the limits m — t- cxo and m — ?• are detailed in Sects. 5.2 and 5.3, respectively, 
while the floating fixed point scenario with finite m is elucidated in Sect. 5.4. 

5.1 Analytic structure of the beta functions 

The key feature for understanding the RG flow implied by (4.40) is the analytic structure of the 
beta functions (4.41). Eq. (A. 14) illustrates that summing up the tower of Matsubara states 
either leads to hyperbolic functions which are well-defined on the entire real axis or trigonometric 
functions which give rise to poles for certain values x. Which situation is actually realized is 
determined by the relative sign between the kinetic ((9^)-terms and the potential terms including 
the spatial momenta and cosmological constant. 

In an Euclidean theory with conventional (healthy) kinetic terms the n? and terms come 
with a relative plus sign, refiecting that the "time" -like and spatial momenta stem from a positive 
semi-definite Laplace operator. In this case the Matsubara sums solely give rise to hyperbolic 
terms. Performing the Wick rotation to Lorentzian signature changes the relative sign and thus 
one encounters trigonometric terms. Essentially, this reflects the fact that in the two settings 
the symbol of the differential operator appearing in the definition of the kinetic term is positive 
definite or of mixed signature. 

For the truncation (4.3) this simple picture does not hold. Inspecting the threshold functions 
entering into (4.41), a careful analysis of the signs of the roots (A. 21) indicates that the analytic 
structure of the beta functions depends on the value of the dimensionless cosmological constant 
A. As indicated in Table 2 there are three distinguished regions in the g-X plane. For A < — 4Co 
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the beta functions for Euclidean and Lorentzian signature are build only from hyperbolic and 
trigonometric functions, respectively. For A > 1/2 the analytic properties are interchanged, i.e., 
the Euclidean signature leads to trigonometric terms while the beta functions originating from 
Lorentzian signature are purely hyperbolic. In the intermediate region one always encounters a 
mix of hyperbolic and trigonometric terms. 

The latter property can be traced back to the structure of the functional traces (4.23) con- 
tributing to the flow. For Euclidean signature, the transverse-traceless fluctuations give rise to 
hyperbolic terms while the conformal mode contributes the trigonometric pieces. Ultimately, 
the occurrence of the trigonometric terms is a consequence of the well-known conformal factor 
problem of gravity:^ Performing the second variation of the Ricci scalar leads to "wrong-sign" 
spatial momenta for the conformal scalar h. Combined with the standard "time" derivatives of 
the ansatz (4.3) this results in a "wrong" relative sign between kinetic and potential contribu- 
tions. In other words the Euclidean conformal sector essentially contributes as a standard scalar 
field in a Lorentzian signature spacetime. As a consequence, the beta functions develop poles at 
values X, where a given Matsubara mode leads to a vanishing denominator in (A. 14). 

5.2 The compactification limit 

We now study the flow of the beta functions (4.40) for /c-independent T in the limit of a collapsing 
time circle, T — t- 0. In this limit, the Matsubara modes become infinitely heavy and decouple, 
so that the flow becomes essentially d-dimensional. In order to account for this effect we trade 
the D-dimensional Newton constant appearing in the Einstein-Hilbert ansatz (4.3) for its 

(i-dimensional analog 

At the level of the dimensionless couplings (4.31), this relation translates into 

mgi''^ =27Tgi'K (5.2) 

The beta functions governing the scale dependence of g^^ can then be deduced by first rewrit- 
ing (4.40) in terms of the d-dimensional Newton constant and subsequently taking the limit 

*The conformal factor problem is also present in the metric formulation of the FRGE. In this framework the 
wrong-sign kinetic term of the h field is compensated by a suitable choice of regulator function TZk . The structure 
of the flow equation (1.1) then ensures that the conformal mode contributes to the flow exactly like a scalar field 
with positive-definite kinetic term. 
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m — )• oo. Utilizing the intermediate results of Appendix B.l we obtain 

a^^f = /3('^)(5^'^\ A) , dtXk = /3a(<7^'\ A) (5.3) 

where 

/3f)=(d-2 + #)5('^), 



/3, = [v7 - 2j A + j^-^ [d,^ + 1 j («>^/y-2A) - ir?V"^$i;2(-2A) 
The anomalous dimension of the d-dimensional Newton constant is given by 

''^ " l + 27r5W 52(A) ^ ^ ^ 

where the functions Bi and B2 have simplified to 

^1 - I 6 (^2T + 1) «d/2_i - MA^d/2 + ( MA " C«2tC2T —)^d/2 



R 4 f 1 / J I nrfsl.O d-2fE,l,0 /d-2 , ^ 2(d-3) \ 2 2,0 
^2 = 7^3d72 6 (^2T + 1) 'l^d/2~l - 2dX^d/2 + ( MA " f^2TC2T 



(5.6) 

(47r)d/2 I 6 V"-^!' -r -^d/2-1 - 2dX-^d/2 V MA " "^T^^l —J ^ d/2 

and all threshold functions are evaluated at = — 2Afc. 

At this stage the following remarks are in order. Inspecting (5.4) we first observe that the part 
of the beta functions induced by the mass dimension of G^^^ matches the one expected for the 
Newton constant in d dimensions. Moreover, the signature dependence of the beta functions has 
dropped out in the compactification limit. This is expected, since all information on the signature 
of spacetime has been encoded in the T circle. Finally we note that the quantum contributions 
to the running of the cosmological constant and the Newton constant are essentially provided by 
the transverse-traceless fluctuations on the d-dimensional spatial slice and one scalar field. This 
is precisely the off-shell-field content of d-dimensional gravity upon gauge fixing diffeomorphism 
invariance. In this light, the decoupling of the transverse spatial vector observed in (4.22) seems 
necessary to reproduce the correct T — ?• limit of the gravitational flow equations. 

A crucial ingredient in understanding the RG flow implied by the system (5.3) are the fixed 
points of the RG flow which by definition satisfy Pg{gi'^\ X*) = 0, Px{9i^\ X*) = 0. The beta 
functions posses a Gaussian fixed point (GFP) gi'^^ = 0, A* = 0. This fixed point corresponds to 
the free or Gaussian theory. Moreover, the beta functions also give rise to a unique non-Gaussian 
fixed point (NGFP) with gi'^^ > 0, A* 7^ for any dimension d > 3 whose position can easily be 
found numerically. For 3 < d < 7 their values together with the universal scaling variable 

_ - / 2/(^^-2) 

n = A* [gi ' j (5.7) 
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d 


(d) 

g* 


A* 




^1,2 


3 


0.24 


0.30 


0.02 


0.89±3.22i 


4 


0.77 


0.28 


0.22 


2.69±4.63i 


5 


3.17 


0.29 


0.62 


4.55 ±6. 26i 


6 


15.3 


0.29 


1.14 


6.64 ± 7.86i 


7 


84.0 


0.30 


1.75 


8.97±9.41i 



Table 3: Position and critical exponents of the non-Gaussian fixed point of the beta functions 
(5.4) for selected spatial dimensions 3 < < 7. The universal product r*, eq. (5.7), and the 
critical exponents of the NGFP are very similar to the ones found in the metric Einstein-Hilbert 
truncation [10, 54]. 

are summarized in Table 3. Given a fixed point of the flow equation, the linearized flow near 
the fixed point is governed by the Jacobi matrix B = (Bij) = dj/3i\^^^^_^(d) . The stability 
properties of the fixed point are then encoded in the stability parameters 9 given by minus the 
eigenvalues of B. For the NGFP of our system these are given in the last column of Table 3. 
The positive real part of 9 thereby indicates that the NGFP is UV-attractive in both gi'^^ and A* 
while the non-zero imaginary part implies that RG trajectories coming close to the fixed point in 
the UV will spiral into the NGFP as /e — )• oo. 

Remarkably, the properties of the NGFP reported in Table 3 are very similar to the NGFP 
featuring in the Asymptotic Safety program for metric gravity. In particular the universal scaling 
variable and the critical exponents turn out to be similar to the typical values obtained in the 
analysis of the beta functions of the Einstein-Hilbert action in the metric formalism [10, 54]. Thus 
it is tempting to speculate that the UV behavior of quantum gravity in the metric formulation 
and the ADM-decomposed formulation at T = are actually governed by the same universality 
class. 

The phase diagram obtained by integrating the beta functions (5.4) numerically is shown in 
Fig. 1. The UV-behavior of all trajectories with positive Newton constant is governed by the 
NGFP. Tracing the RG flow towards the IR, there are three classes of trajectories: the separatrix 
(Type Ila) connects the NGFP with the GFP. Trajectories which pass to its left (right) are 
classified as Type la (Type Ilia), respectively, and lead to a negative (positive) IR- value of the 
cosmological constant. The flow pattern is virtually identical to the one obtained in the metric 
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-0.4 -0.2 0.2 0.4 -0.4 -0.2 0.2 0.4 



Figure 1: Phase diagram originating from the foliated Einstein-Hilbert truncation in the Hmit 
T — )• in (i = 3 (left) and d = 4 (right). Along the red line the anomalous dimension (5.5) 
diverges. The flow is strikingly similar to the one obtained within the metric Einstein-Hilbert 
truncation [53]. 

Einstein-Hilbert truncation [53]. 

5.3 The decompactification limit 

In the last subsection, we found that in the limit T — )• all signature effects drop out of the 
beta functions (4.41). In this subsection we analyze the opposite limit where the foliation be- 
comes infinitely extended T — ?• oo. In this limit, the discrete Matsubara sums in (A. 14) become 
continuous integrals. As a consequence, the trigonometric terms in the beta functions develop 
a branch-cut singularity leading to diverging expressions. Thus the beta functions make sense 
in the case where they consist of hyperbolic terms only. Inspecting Table 2 shows that the beta 
functions are well-defined in the regions A < — 4Co in Euclidean and A > 1/2 for Lorentzian 
signature only. In particular they become ill-defined in the central region — 4Co < A < 1/2 for 
both signatures, since there they always include both hyperbolic and trigonometric terms. 

From Table 3 we note that the NGFP in the compactification limit is situated in the central 
region. As a consequence we do not expect to find a NGFP for the T — ?• oo limit. Indeed a 
numerical analysis confirms that there is no NGFP for positive Newton constant. This should, 
however, not be interpreted as a failure of Asymptotic Safety. As discussed in Subsection 5.1 
the appearance of the trigonometric terms and thus the break down of the beta functions in the 
central region originates from the conformal-factor problem intrinsic to the ansatz (4.3). Clearly, 
our truncation is too simple to allow for a dynamical resolution of this problem. Addressing this 
problem systematically presumably requires the inclusion of higher derivative operators like the 
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square of the intrinsic Ricci scalar {^'^^R)'^ in the truncation ansatz. Since the inclusion of such 
higher-derivative terms in the truncation ansatz is known to be very laborious [55, 9, 56] and 
beyond a first investigation of foliated RG flows, we do not pursue this direction further and 
leave it to future research. 

5.4 RG flows at finite m 

In our final analysis we complement the system (4.40) by a beta function for the Matsubara mass 
mfc. Our basic assumption is that the dimensionless also approaches a finite value at the 
NGFP. The supplementary beta function, capturing this behavior is 

dtnik = (3m{9A,m) (5.8) 

where the NGFP is a common zero of all three beta functions 

/3g{g*,K,m^) = , l3x{g*,X*,'m^) = , f^mig*, K,m*) = (5.9) 

with m* 7^ 0. Since the computation of f3m is beyond the scope of the present work, we will 
approximate r/ifc = 7^ as a free parameter and investigate the parametric dependence of 
(4.40) on m = m*. 

An immediate consequence of approximating m = const, dtrrik = is that the dimensionfull 
length of the foliation direction scales with k according to its canonical dimension T (x k~^. In 
[35] this was dubbed the "floating fixed point scenario" . In this scenario the time-interval length 
is inverse proportional to the RG scale k. This means that the time circle decompactifies in the 
IR as one would expect from a physically viable theory. On the other hand in the deep UV the 
circle collapses in a controlled way, so that the presence of the "extra dimension" and signature 
still enters into the beta function through the dimensionless parameter m. 

The possibility of this additional k dependence has already been anticipated in the definition of 
the anomalous dimension of the Newton constant, (4.32), so that the beta functions (4.41) remain 
valid in this case as well. For m > the beta function for both signatures remain regular 

throughout the entire central region. When giving numerical results we will choose m = 27r for 
illustration purposes. 

We analyze the m dependence of the NGFP arising from the system (4.41) in d = 3, corre- 
sponding to a Z) = 3 -|- 1-dimensional spacetime. In Subsection 5.2 we established the existence 
of the NGFP in the limit m — )■ 00. In this case the Euclidean and Lorentzian NGFP fall on top 
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Figure 2: Parametric m dependence of the position {g^,, A*} of the NGFP (upper left panel) and 
its complex stability coefficients 612 = 0' i: iO" (upper right panel). The corresponding location 
in terms of the d-dimensional Newton constant (5.2) is shown in the lower panel. The blue 
circles and magenta squares correspond to Euclidean and Lorentzian signature respectively. In 
the compactification limit m — t- 00 the Lorentzian and Euclidean NGFP coincide. 



of each other. In terms of the dimensionless Newton constant in D dimensions g\ = g* it is 
situated at {gi^\ X*} = {0,0.3}. The parametric m dependence of the NGFP is shown in Fig. 
2. Lowering m, the signature effects drive the Euclidean and Lorentzian NGFP apart. However, 
the signature dependence is rather mild, so that its location and critical exponents are rather 
robust with respect to a change in m. In particular we observe that the real part of the critical 
exponents is always positive so that the NGFP is UV attractive in both (7* and A* for all values 
m. Thus the NGFP constitutes a viable candidate fixed point for the gravitational Asymptotic 
Safety program in both Euclidean and Lorentzian signature spacetimes. 

We illustrate these findings by giving the numerical results for the position and critical ex- 
ponents of the NGFP at m = 27r in Table 4. In this case the signature affects the value of the 
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e 


9* 


A* 


5* A* 


^1,2 


+1 


0.23 


0.30 


0.070 


0.92 ± 3.27i 


-1 


0.24 


0.30 


0.072 


0.87±3.16i 



Table 4: Fixed point values and critical exponents in the floating fixed point scenario for m = 2tt 
and D = 3 + l. 

universal product 5f*A=i, and the real and imaginary part of the critical exponents at the 5% level. 
This highlights the remarkable robustness of the flow implied by (4.40) with respect to changing 
the spacetime signature. 

This rather remarkable stability also extends away from the NGFP. This feature is illustrated 
in Fig. 3 where we constructed the phase diagrams for Euclidean and Lorentzian signature flows 
by numerically integrating the flow equations for the D-dimensional Newton constant and 
the cosmological constant A^ for m = 2tt. In both cases the flow is dominated by the interplay 
between the NGFP in the UV and the Gaussian fixed point located at the origin {g^,, A*} = {0, 0} 
in the infrared. The flow diagrams turn out virtually identical, consolidating the assessment that 
the signature of spacetime plays a subleading role when determining the gravitational RG flow. 
Moreover, they coincide with their counterparts obtained in the compactification limit displayed 
in Fig. 1. Furthermore the corresponding classification of RG trajectories is the same as the one 
obtained from the flow equations in the metric formulation [53] . 

6 Conclusions and outlook 

In this work we constructed the Wetterich equation [24] which captures the non-perturbative 
RG flow of the effective average action F^ on the theory space of projective Hofava-Lifshitz 
gravity [18]. Its solutions interpolate continuously between the classical action S = T^^oo and 
the full quantum effective action F = T^^q. In contrast to the metric construction [7], this 
novel functional renormalization group equation (FRGE) encodes the gravitational degrees of 
freedom in the lapse function, the shift vector and the spatial metric originating from the ADM 
decomposition [26]. The construction imprints a foliation structure on the quantum spacetime, 
giving rise to a Euclidean "time" direction. Owed to the foliation, the flow equation closely 
resembles the ones encoding the RG flows of quantum field theory at finite temperature [33]. 
Analogously to the metric construction [7], our FRGE for the ADM fields employs the back- 



32 




A 



Figure 3: Flow diagram for Euclidean (left panel) and Lorentzian (right panel) signature for 
m = 27r in = 3 + 1 spacetime dimensions. 

ground field method with a linear background split and constructs the IR regulator from the back- 
ground quantities. The resulting RG flows explicitly preserve background- foliation-preserving 
symmetries Diff(A^,S) C Diff(A^). This invariance ensures that the RG flow does not gener- 
ate interactions that violate this symmetry dynamically. Thus our equation constitutes a natural 
starting point for studying the RG flow of anisotropic gravity theories like Hof ava-Lifshitz gravity, 
possibly shedding light on hitherto open questions concerning the perturbative renormalizability 
and a possible dynamical restoration of Lorentz symmetry at low energies of the model. 

As explained in Sect. 3 there is a clear relation between the RG flows captured by the metric 
formulation of gravity and projective Hofava-Lifshitz gravity. While, upon gauge-fixing, the off- 
shell field content of the two theories looks identical, the weaker symmetry requirements of the 
anisotropic case entail that the theory space of Hofava-Lifshitz gravity contains the one of metric 
gravity as a subspace. The result obtained in [7], showing that metric formulation of the FRGE 
does not generate non-diffeomorphism invariant interactions dynamically, indicates that the latter 
is closed under RG transformations. 

As a first application of our novel FRGE, we studied the RG flow captured by the ADM- 
decomposed Einstein-Hilbert action (4.3). The beta functions encode the scale dependence of the 
Newton constant and the cosmological constant and parametrically depend on the dimensionless 
Matsubara mass m, associated with the "size" of the "time direction". As a key result, already 
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anticipated in [35], we establish that the RG flow resulting from the foliated flow equation gives 
rise to a non-Gaussian fixed point intrinsic to the gravitational Asymptotic Safety scenario. This 
result is independent of the signature of spacetime, as long as m is sufficiently large. Moreover, 
the phase portraits shown in Figs. 1 and 3 obtained from the foliated flow equation are strikingly 
similar to the ones found within the metric Einstein-Hilbert truncation [53]. Based on this result, 
we expect that the foliation structure of spacetime should not play a crucial role when determining 
the short distance behavior of the theory. 

Naturally, there are a whole range of applications that can be explored based on the FRGE 
(3.34). Obviously, expressing the gravitational degrees of freedom in terms of the ADM fields 
together with the invariance under background-foliation-preserving diffeomorphisms prepares the 
stage for studying the RG flow of Hofava-Lifshitz gravity. The truncation ansatz (4.3) can 
straightforwardly by supplemented by interaction monomials which break DifF(7W) but are in- 
variant under Diff(A^,S). Interesting examples could include a higher-derivative term in the 
intrinsic curvature / (Pxdt ^/aN \^'^^R^ and the square of the Cotton tensor J d^xdt ^/o^NCijC^^ 
which are at the heart of the perturbative renormalizability of the theory. Moreover, the inclusion 
of matter fields with anisotropic couplings is straightforward. The simplest extension which in- 
cludes a scale-dependent coupling in the extrinsic curvature terms is currently under investigation 
and the results will be reported elsewhere [36]. 

The presence of the foliation structure also constitutes an important step toward connect- 
ing the continuum FRGE to the Monte-Carlo simulations of the gravitational partition function 
within the Causal Dynamical Triangulations (CDT) program [31]. In a sense our novel flow equa- 
tion contains the continuum analogue of the foliation structure that is supposed to be responsible 
for the well-defined classical limit in CDT. Following the ideas [57, 58, 59, 56] one could, e.g., use 
the foliated flow equation to study diffusion processes on spatial slices and examine if these can 
be matched to the CDT data recently reported in [60]. This may shed new light on the question 
whether the continuum limit of CDT corresponds to an isotropic or anisotropic gravity theory as 
suggested in [61, 62]. 
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A Trace technology 

In this appendix, we collect the technology and intermediate results used for evaluating the 
functional traces (4.24) and (4.25). We start with the collection of some standard heat-kernel 
techniques for evaluating the operator traces on the spatial slices in subsection A.l while the 
Matsubara sums are discussed in subsection A. 2. 

A.l Traces on the spatial slice 

The right hand side of the flow equation (4.23) naturally gives rise to operator traces tr[VF(A)], 
restricted to the spatial slices and containing functions of the spatial Laplacian A = —a^^DiDj. 
Introducing the Laplace anti-transform W{s) 



oo 



W{z) = / dsW{s)e-^' (A.l) 

JO 

these operator traces can be written as 

roc 

tr[VF(A)]=/ dsW {s) tr [e-'^] . (A.2) 
Jo 

The trace of the heat kernel has the well-known asymptotic expansion for s — )• 

tr [e-^] = ^^-i^ j d'xV^ [d, + lsRd. + .. .] (A.3) 

where di is shorthand for the dimension of the internal space and the dots denote curvature 
monomials build from higher orders of the spatial curvature tensors. 

Substituting the expansion (A.3) into (A.2) leads to a curvature expansion of the operator 
traces in terms of the intrinsic curvatures 



tr [W{A)] = J d'^xV^ ds W{s) [d, + \sRd,+ 

d'^xV^ [d, Qd/2[W] + \d^R Qd/2-i m + . 



(47r)^/2 

with the Q-functionals defined as 



(A.4) 



/•oo 

Qn[W] = / dss-'^Wis). (A.5) 
Jo 
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For n > 0, the definition of tlie F-functions allows to re-express these functionals as a Mellin 
transform 



1 



oo 



Qn[W] = -^l dzz^"^W{z). (A.6) 

In order to write down the beta functions, it is useful to rewrite these Q-functionals in terms 
of dimensionless threshold functions. The VF-functions appearing in the traces (4.24) and (4.25) 
are of the form 

This functional form motivates defining 

1 p , i?(0)(z)-zi?(0)'(z 



^""^ - Y{n) j, {z + Ri^){z) + wY [z + Ri^){z)y' 

r dzz--^ ^^'^^'^ 



(A.8) 



{z + R'S>){z) + wY {z + Ri^){z)) 



satisfying the recursion relations 



A$P-'?(^) = -p^p+^^'i{w) , ^4rH = -P^r^'^M . (A.9) 

The definitions (A.8) then imply 

Q„ [WP''^] = 2fe2(n-p-9+l) $P.'?(u;/fc2) , Q„ Iw^'l] = fe2("-P-9+l) $P.9(-«;/A;2) . (A.IO) 

Inspecting the traces (4.24) and (4.25) one furthermore observes that the functions W and W 
always appear in certain combinations. In this view it is also natural to introduce 



(w) ^ - IrnK'^H (A.ii) 

where t^tv is the anomalous dimension of the Newton constant (4.32). 

Up to this point the profile function R^^\z) entering the IR cutoff was left unspecified. For 
explicit computations it turns out to be convenient to work with the optimized cutoff [63] 

r(^\z) = {1 - z) e {1 - z) . (A.12) 

For this choice the integrals appearing in the threshold functions (A.8) can be carried out ana- 
lytically. The result degenerates in q and reads 

We will implement this choice of cutoff, when resumming the Matsubara modes in the next 
subsection. 
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A. 2 Resumming the Matsubara frequencies 

Owed to the foliation structure, our flow equation naturally contains infinite sums over the Mat- 
subara modes. For special choices of the regulator, as, e.g., the optimized cutoff (A. 12), these 
sums can be carried out analytically. For the threshold functions appearing in the main text, this 
can be done via the master formula 

oo oo ^ 

^-^ rfl + x"^ X tanhfvrx) ' ^—^ v? — x^ x tanfvrx) 

n=— oo ^ ' n=— oo ^ ' 

related by analytic continuation. Based on this relation the sums in (4.33) and (4.35) can be car- 
ried out explicitly. Note that throughout this subsection we work exclusively with the optimized 
cutoff, such that the threshold functions take the form (A. 13). 
We start with the transverse-traceless sector, by defining 

oo oo 

ih^^n^y ^rM. tr= {h^^r?Y (a.is) 

n=— oo n=— oo 

with W2'Y given in (4.34). Inspecting (4.33), we observe that the transverse-traceless sector con- 
tains functions T with r = only. For p = 1 the sum is of the form (A.14) and yields 



= . ' (A.16) 



r(Z + 1) m V 1 - 2A t^j^h (^v^e(l-2A) 



The T functions with higher values of p are found by applying the recursion relation (A. 9). 
Exploiting that T^'*^ = \d\ T^^'*^, we obtain 

/ 



^2,0 _ 1 s/l-n 

I ~ T(l+1) 2m(l-2A) 



V 



, ,^ ^ + , I . (A.17) 

VI - 2A tanh (^Ve(l-2A)) m sinh (^Ve(l-2A)) 

Resumming the threshold functions in the scalar sector is slightly more involved, but can by 
carried out along the same lines. Analogously to (A. 15), we start by introducing 

oo oo 

{hmWy ^f\wo), ^'f^ (7"^V)'' ^^^(^wo) . (A.18) 

n=— oo 71=— oo 

The wq, given in (4.36), leads to denominators which are (powers of) fourth order polynomials in 
n. This case can be reduced to the master formula (A.14) by first factorizing the denominators 
into second order polynomials and subsequently carrying out a partial fraction decomposition. For 
the purpose of computing the ^-functions recursively, it is useful to generalize the denominator 
1 + Wo by introducing auxiliary parameters /3, 7 

l + u;o = 7-2A + /3^imV-^imV . (A.19) 

/3=7=1 
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The polynomial in n factorizes into 

l + ^o = -(^ 

with 



_ HL 
2 A 



w± 



^Aw^ i[w^) +^a(^-'^) 



(A.20) 



(A.21) 



=7=1 



Substituting (A.20) into ^'f^ and carrying out the partial fraction decomposition casts the sum 
into the form (A. 14), so that it is easily evaluated to be 

1 d-2 ttX / 1 1 



1,0 



(A.22) 



' r(/ + 1) d — 1 {w^ — w_) \^w^ tan{TT^/ui^) yJwZ tan(7r^t(;_) ^ 
Based on this expression the other ^'^'^ can be found recursively by taking derivatives with 



respect to the auxiliary parameters introduced in (A. 19). In particular the ^^f^ satisfy 



(A.23) 



=7=1 



For p = 1 this yields 



2,0 



1 



d-2 



tA2 



"r(/+l) 1 J (w+—w-.)'^ '] W+—W- 

1 f 1 



1 



1 



+ 



2ui+ y^/w^ ta.n{TT^/W+) ~^ sin'' {n ^^^U^ J ' 2u)- 

Furthermore, powers of in the numerator can be generated by taking derivatives with respect 
to /3. For the ^'^'^ appearing in (4.38) this implies 

-.2,1 



w+ tan(7r^?i)+) ^ui_ tan(7r^«)- ) 
1 f 1 



+ 



w- tan(7r^ui_) ~'~ sii?(7r77^3) 



(A.24) 



^^f'^ = 

' 3d - 4 5/? ' 



(A.25) 



/3=7=1 



which evaluates to 

.2,1 1 



{d - 2)2 ttA + 
r(/ + 1) 2((i - l)(3d - 4) ^ (u7+ - W-Y 

1 / 3ui++to_ 3ui_+ui+ 



The corresponding expressions for the sums T and ^ are easily obtained from the results given 
above by taking into account the different prefactors (A. 13), so that we refrain from giving explicit 
formulas. 

For the purpose of writing the beta functions in a compact way, it is furthermore useful to 
introduce the analog of the q functionals (A. 11) setting 



+ 



(A.26) 



rpP,r _ 



p,r 1 



p,r 



(A.27) 
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where the T and 5 are reminders that the threshold functions appear in the transverse-traceless 
and scalar traces, respectively. Given these results the evaluation of the functional traces (4.24) 
and (4.25) is now straightforward. 



B Limits of the threshold functions 

In order to understand the structure of the flow implied by the beta functions (4.40), it is useful to 
study the two limiting cases where the time circle collapses, T — t- 0, or becomes infinitely extended, 
T — )■ oo. In Sect. 5 these two cases are called the compactification and decompactification limit 
respectively. This appendix then collects the corresponding limits of the threshold functions 
constructed in Sect. A. 2. 

B.l The compactification limit T — )• 

In the limit T — t- 0, where the time circle collapses, the Matsubara mass (4.31) becomes infinite, 

lim m = oo . (B.l) 

As a consequence the massive Matsubara modes decouple from the flow. The only term contribut- 
ing to the Matsubara sums (A. 15) and (A. 18) is the zero mode n = 0. It is then straightforward 
to verify that 




limTf''^= { °' (B.2) 

r/0<p-l 



and 



hm ^f^' 



$f(-2A), r = 0, 

^ (B.3) 
, r / < 2p - 1 . 



Taking the compactification limit thereby commutes with carrying out the Matsubara sums. In 
particular taking the limit m — )• oo in eqs. (A. 16), (A. 17), (A. 22), (A. 24) and (A. 26) confirms 
the general identities (B.2) and (B.3). This provides a non-trivial crosscheck for the threshold 
functions T and ^. As in the previous section, the result (B.2) and (B.3) also holds for T and ^ 
if the r.h.s. is replaced by $^(— 2A). 
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B.2 The decompactification limit T — oo 

The limit T — )• oo is slightly more involved. In this case the Matsubara mass vanishes 

lim m = (B.4) 

T— >oo 

SO that the discrete sums essentially turn into a continuous integral. Inspecting the threshold 
functions one observes that this limit gives rise to finite results if and only if the threshold 
functions are of hyperbolic nature. Table 2 then indicates that this is the case in two regions on 
the g-X plane 

I : A < -4Co , e = 1 

(B.5) 

II: A>l/2, e = -l. 

Outside these regions the beta functions oscillate infinitely rapidly and are thus ill-defined. There- 
fore we will restrict ourselves to giving the limit behavior of the threshold functions appearing in 
the beta functions (4.40) for these two regions only. To simplify notation, it then turns out to be 
useful to extract the prefactors from (A. 21) and define 



. /3(3d-4) , (m^\\d-2l 

''^=W^ nW^^ +:^t(7-2A) 



d-lX 



(B.6) 

/3=7=1 



which are manifestly positive in both regions. 

Taking the limit m — )• in the tensor sector, (A. 16) and (A. 17), is straightforward and yields 

i^^o = ^riri) (1 - 2A)p-i \r:w^ ' ^ = i'2- (^-^^ 

The scalar sector is slightly more involved. Here we first note that since e = 1, A < in case I 
and e = — l,A>Oin case II, the arguments of the square-roots appearing in (A. 22), (A. 24) and 
(A. 26) are actually negative which requires the analytic continuation of the formulas to hyperbolic 
functions. Subsequently taking the limit (B.4) in (A. 22) gives 



lim m = --1^ '-^ ' ' - ^ • (B.8) 

The higher order threshold functions can then again be obtained via the recursion relations (A. 23) 
and (A. 25). In particular 



lim mM'2'°=^ i_ i^^/__4_ (i ^] + ^ + ^ 

m^O ' r(Z+l) \d-l) 2|A|3/2 {W+-W.Y | w+-w- \^/^ ' ^ ^ 



+ 



urn m^^ - (u,+ -w.Y \ 
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(B.9) 



Again the corresponding expressions for the sums T and ^ are easily obtained from the results 
given above by taking into account the different prefactors (A. 13), so we refrain from giving 
explicit expressions. 
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